home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / MINFIT.z / MINFIT
Text File  |  1996-03-14  |  4KB  |  133 lines
  1.  
  2.  
  3.  
  4. ____MMMMIIIINNNNFFFFIIIITTTT((((3333FFFF))))                                                        ____MMMMIIIINNNNFFFFIIIITTTT((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      MINFIT, SMINFIT  -  EISPACK routine.  This subroutine determines, towards
  10.      the solution of the linear
  11.         T system AX=B, the singular value decomposition A=USV  of a real
  12.         T M by N rectangular matrix, forming U B rather than U.  Householder
  13.      bidiagonalization and a variant of the QR algorithm are used.
  14.  
  15.  
  16. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  17.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee  mmmmiiiinnnnffffiiiitttt((((nnnnmmmm,,,, mmmm,,,, nnnn,,,, aaaa,,,, wwww,,,, iiiipppp,,,, bbbb,,,, iiiieeeerrrrrrrr,,,, rrrrvvvv1111))))
  18.           iiiinnnntttteeeeggggeeeerrrr          nnnnmmmm,,,, mmmm,,,, nnnn,,,, iiiipppp,,,, iiiieeeerrrrrrrr
  19.           ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn aaaa((((nnnnmmmm,,,,nnnn)))),,,, wwww((((nnnn)))),,,, bbbb((((nnnnmmmm,,,,iiiipppp)))),,,, rrrrvvvv1111((((nnnn))))
  20.  
  21.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssmmmmiiiinnnnffffiiiitttt((((nnnnmmmm,,,, mmmm,,,, nnnn,,,, aaaa,,,, wwww,,,, iiiipppp,,,, bbbb,,,, iiiieeeerrrrrrrr,,,, rrrrvvvv1111))))
  22.           iiiinnnntttteeeeggggeeeerrrr          nnnnmmmm,,,, mmmm,,,, nnnn,,,, iiiipppp,,,, iiiieeeerrrrrrrr
  23.           rrrreeeeaaaallll             aaaa((((nnnnmmmm,,,,nnnn)))),,,, wwww((((nnnn)))),,,, bbbb((((nnnnmmmm,,,,iiiipppp)))),,,, rrrrvvvv1111((((nnnn))))
  24.  
  25.  
  26.  
  27. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  28.      On INPUT
  29.  
  30.      NNNNMMMM must be set to the row dimension of two-dimensional array parameters
  31.      as declared in the calling program dimension statement.  Note that NM
  32.      must be at least as large as the maximum of M and N.
  33.  
  34.      MMMM is the number of rows of A and B.
  35.  
  36.      NNNN is the number of columns of A and the order of V.
  37.  
  38.      AAAA contains the rectangular coefficient matrix of the system.
  39.  
  40.      IIIIPPPP is the number of columns of B.  IP can be zero.
  41.  
  42.      BBBB contains the constant column matrix of the system if IP is not zero.
  43.      Otherwise B is not referenced.  On OUTPUT
  44.  
  45.      AAAA has been overwritten by the matrix V (orthogonal) of the decomposition
  46.      in its first N rows and columns.  If an error exit is made, the columns
  47.      of V corresponding to indices of correct singular values should be
  48.      correct.
  49.  
  50.      WWWW contains the N (non-negative) singular values of A (the diagonal
  51.      elements of S).  They are UNORDERED.  If an error exit is made, the
  52.      singular values should be correct for indices IERR+1,IERR+2,...,N.
  53.         T
  54.  
  55.      BBBB has been overwritten by U B.  If an error exit is made,
  56.         T the rows of U B corresponding to indices of correct singular values
  57.      should be correct.
  58.  
  59.      IIIIEEEERRRRRRRR is set to Zero       for normal return, K          if the K-th
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ____MMMMIIIINNNNFFFFIIIITTTT((((3333FFFF))))                                                        ____MMMMIIIINNNNFFFFIIIITTTT((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      singular value has not been
  75.         determined after 30 iterations.
  76.  
  77.      RRRRVVVV1111 is a temporary storage array. Calls PYTHAG(A,B) for sqrt(A**2 +
  78.      B**2).  Questions and comments should be directed to B. S. Garbow,
  79.      APPLIED MATHEMATICS DIVISION, ARGONNE NATIONAL LABORATORY
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.